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Abstract 



We consider time-changed Poisson processes, and derive the governing difference-differential 
equations (DDE) these processes. In particular, we consider the time-changed Poisson pro- 
|V | . cesses where the the time-change is inverse Gaussian, or its hitting time process, and discuss 

^ ■ the governing DDK's. The stable subordinator, inverse stable subordinator and their iterated 

versions are also considered as time-changes. DDE's corresponding to probability mass func- 
tions of these time-changed processes are obtained. Finally, we obtain a new governing partial 
differential equation for the tempered stable subordinator of index < /3 < 1, when /3 is a 
rational number. We then use this result to obtain the governing DDE for the mass function 
. of Poisson process time-changed by tempered stable subordinator. Our results extend and 

in 
p 

ly-^ I Key words: Hitting times; inverse Gaussian process; stable processes; time-changed process. 



complement the results in Baeumer et al. [3] and Beghin et al. [1] in several directions. 

Key words: Hitting times; inverse Gaussian process; stable processes; tin 
subordination; tempered stable processes; difference-differential equation. 



; 1 Introduction 



Recently there has been an increasing interest to consider time-changed stochastic processes 
that yield solutions of fractional Cauchy problems, or solutions of higher order partial differen- 
tial equations (PDE). The factional Cauchy problems can be used to model various phenomena 
in a wide range of scientific areas including physics, telecommunications, turbulence, image 
processing, biology, bioengineering, hydrology and finance, see [Sj HI [TTl [181 [19] . There is an in- 
teresting connection between continuous time random walks and fractional Cauchy problems, 
see [201 [HI [19] 

It is well known that the Poisson process N{t) with parameter A > solves the following 
difference- differential equation (DDE) 

^PfcW = -Apfe(t) + Apfc_i(t), (1.1) 
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where Pk{t) = F{N{t) = k} is the probabihty mass function (pmf). The fractional Poisson 
process A^^(t), which is a generahzation of the Poisson process, solves the following fractional 
DDE (see e.g. Laskin (2003)) 

^Pi{t) = -Xpi{t) + Xpl,{t), (1.2) 

where (t) = F{N^{t) = k} and 

d^u{t,x) 1 du{r,x) dr 



;i.3) 



dt^ r(l - (3) Jq dr {t - ry 

for < /3 < 1, denotes the Caputo fractional derivative (see e.g. Caputo (1967)). Let 
u{s,x) = e~'^^u{t,x)dt = Ct{u(t,x)) be the Laplace transform (LT) of u(t,x) with respect 
to the variable t. Then Laplace transform is given by 

A (^^0^) = ^'^(^' ^) - ^'"'«(0' (1-4) 

The inverse Gaussian (IG) process G{t) (or IG subordinator) has been found useful in 
financial modeling and is defined by (see Applebaum, 2009, p. 54) 

G{t) =mi{s > 0;B{s) +-fs > 5t}, (1.5) 

where B{t) is the standard Brownian motion. Note that G{t) ~ lG{6t, 7), the inverse Gaussian 
distribution with density 

(7(x,t) = (27r)-^/'(5t)x~=^/2g57t-|(4^+7^x)^ 3;>0. (1.6) 

Also, note that 

H{t) = inf{s > : G{s) > t} (1.7) 

denotes its first hitting time process. 

In this paper, we consider the Poisson process time-changed by G{t) and by the process 
H{t), and investigate their properties. We consider also the problem of time-change by stable 
and tempered stable process, and derive the underlying PDE's. We only treat the case where 
the process and the time changes are assumed to be independent. Our results extend and 
complement the results in Baeumer et al. [3] and Orsingher et al. [1] in several directions. 



2 IG and its hitting time process as time-changes 

Let N(t) be the Poisson process and G(t) be the IG subordinator. First we consider the 
time-changed process N{G{t)). The probability mass function pk{t) = F{N{G{t)) = k) of the 
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time-changed process N{G{t)) is obtained by the standard conditioning argument as 

Jo k\ 

where Ki,[z) is the modified Bcssel function of third kind with index v, defined by (see e.g. 
Abramowitz and Stegun (1992)) 



1 r'^ 

2 Jo 



(2.2) 



Since, as t ^ oo, a.s. and ^, a.s. (see e.g. Bertoin (1996), p. 92), we have 

t-^oo t X-f ^ ^ 

or equivalently N(G{t)) ~ i^t, a.s. This shows that the subordinated process does not ex- 
plode in any finite interval of time. 

Using the result (see J0rgenson (1992)) 



V TT 7 



we easily obtain 



ENiG{t)) = — and Var(Ar(G'(t))) = — + X\[I{St){-f/'e'^'Ky2{S^t) - {—)\ 
7 7 V TT 7 7 

Remcirk 2.1. For t/ie particular case S — 1 and j — 0, 

F(N(G(t)) = k)= r ^—^^-L=te-''/^-dx 
^ ^ ^ Jo k\ V2^ 

Jo k\ 



-fY{y)dy, 



k\ 

d 



where fviy) = —^=e ^I'^y ^ y > 0. Hence N(G(t)) = N(Xt'^Y), a mixed Poisson process 
evaluated at time t^ . 
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Remark 2.2. Note that N{G{t)) is not a renewal process, since here G{t) is not an inverse 
suhordinator (see Kingman (1964); Grandell (1976)). Indeed, it is the hitting time of B{s)+js, 
which is not a suhordinator. 

The density function g{x,t) of G{t) solves (see Kumar et al. (2011)) 

^g{x, t) - 2S^§^g{x, t) = 26'^g{x, t). (2.4) 
We have the following result. 

Proposition 2.1. The pmf pk{t) of the subordinated process N{G{t)) solves the following 
DDE: 

— Pfc(t) - 257-pfe(i) = 25^X\pk{t) -Pk-iit)]. (2.5) 

Proof. We have 

poo 

Pk{t) = pk{x)g{x,t)dx. (2.6) 

^0 

This implies by dominated convergence theorem 

^Pk{t)^J^ Pk{x)^g{x,t)dx and ^pk{t) ^ Pk{x)^g{x,t)dx. (2.7) 
Using the fact limj;_>.oo 5'(a;, t) = = \im.x^Qg[x.it).i we obtain 



roo Q 

= 25^ / pk{x)—g{x,t)dx 
f°° d 

= —pk{x)g{x,t)dx 

POO 

= -25^ / [-Xpk{x) + Xpk-i{x)]g{x,t)dx 
Jo 

^25^X\pk{t)-pk-i(t)], 



proving the result. □ 

Next we consider the subordination of hitting time of the process G{t). The right continuous 
first hitting time of the process G{t) is defined by 

H{t) = mi{s > : G{s) > t}. (2.8) 

The process H{t) has monotonically increasing continuous sample paths and it is not a Levy 
process. The process N{H{t)) is a renewal process and the inter arrival times follows a 
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tempered Mittag-Leffler distribution. Using Theorem 4.1 of Meerschaert et al. (2010), the 
process N{H{t)) is a renewal process whose iid waiting times J„ satisfy 

F{Jn >x) = ^(e-^^(^)), (2.9) 

and with LT 

E(e-^'") = , r^-— -• (2.10) 

A + 5(^7^ + 2s -7) 

For 6 = 1/ \/2, the rhs of the above equation reduce to 

A 



X + (s + ay -aP' 

where /3 = 1/2 and a = 7^/2. Also, the density of Jn is 



(2.11) 



/,„(x) = ^(x)e— (2.12) 

with 

g(x) = -^[1 - E^i-rixf^)] and 77 = A-a^ (2.13) 
(see Example 5.7 of Meerschaert et al. (2011)). 

Remark 2.3. The condition (5.3) given in Meerschaert et. al.(2010) is satisfied by the process 
G{t), since 



i x\lnx\TT{x)dx = / a;|/na;| 
Jo 



^ e-^'^/'dx 



\/2'KX^ 



< I y-"^\lny\dy 



2t^ JO 

ze ' dz = — = < 00. 



^27r Jo V 27r 

Thus their Theorem 5.2 is applicable on IG density. 

Proposition 2.2. The pmfpk{t) = F{N{H{t)) = k)satisfies 

j^Pk{t) = ^ [(A^ - 25^\)Pu{t) + (-2A2 + 25^\)Pu-i{t) + \^Pk-2{t)\ - 5o(t)Pfc(0) 

+ ^^(0'^K(0), (2.14) 
where ]?fc(0) = —A, A and for /c = 0, 1 and k > 2 respectively. 
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Proof. The density function h{x, t) of H{t) satisfies the following PDE (Vellaisamy and Kumar 
(2011)) 

h{x, t) - 2S-f—h{x, t) = 2S'^—h{x, t) + 25^h{x, 0)5o(t). (2.15) 



dx"^ 

We have for /c > 



dx 



dt 



^Pfe(^) = ^ Pkix)^h{x,t)dx 



Pk[X) 



92 9 
— /i(x,t) -257— /i(x,t) 



dx — 6o{t) I pk{x)h{x,0)dx 




Now 



92 a 
Pkix)—h{x,t)dx = pk{x)—h{x,t) 

d 



00 ^ 
, P'kix)^h{x,t)dx 

00 



Pk{x)—h{x,t) - h{x,t)p'f.{x) 
ox 







(x)/i(x,t)da; (2.16) 



-26-fPk{0)h{0,t) + h{0,t)p'i^Ox) + / p^^^(x)/i(x,t)dx 



since lima;_!.oo /ix(a;, t) = = Mmx^oo h{x,t) (see Vellaisamy and Kumar (2011)). Also 

Pfc(a;)7^/i(x,t)(ix = pfc(x)/i(a;,t)|g° - / p'f,{x)h{x,t)dx 
ox Jo 



-Pk{0)h{0,t) - / p[{x)h{x,t)dx 



(2.17) 



Using the fact hj.{0,t) = 2^7/1(0, t), we get 



d_ 

dt 



-Mi) = ^ 



00 POO 

pf\x)h{x,t)dx + 2S-f / p'^{x)h{x,t)dx + h{0,t)p'^{0) 
io 

-Pfc(0)<5o(t). 



(2.18) 



The result now follows by substituting p^^\x) = X'^[pk{x) — 2pk_i{x) +^^-2(2;)] in fl2.18p . □ 



3 Stable and inverse stable subordinators as time changes 

Let f{x,t) be the density of a /3-stable subordinator D(t) with index < (3 < 1 . Then the 
Laplace transform of D{t) is given by 



E(e-^^(*))=/ e-'^f{x,t)dx = e- 



where is called the Laplace exponent. The density /(x, 1) of D{1) is infinitely differentiable 
on (0, oo), with the asymptotics as follows: (see Uchaikin and Zolotarev (1999)) 

, fl. 2-/9 



/3 



2(1-/3) os,x,- 

f(:^^l)^^ d e-^^-P)^-,) as x ^ 0; (3.1) 

v/27r/3(l - /3) 

/(^'l)~f(rT^j^' asx^oo. (3.2) 
Note that from (13. ip and (13.21) . we have 

lim /(x, 1) = /(0, 1) = anrf lim /(x, 1) = /(cx), 1) = 0. (3.3) 

Let D*{t) = inf{s > : B(s) > t/\/2} denote the first time when the Brownian motion (with 
variance 2t) hits the barrier t. The governing equation corresponding to the density function 
g{x,t) of D*(t) can be obtained by putting 5 = l/\/2 and 7 = in the equation (12.40 . In 
this case, D*(t) is a stable subordinator of index | with Laplace exponent s^/^. We have the 
following proposition. 

Proposition 3.1. Let D2, ■ ■ ■ , -D* be independent 1/2-stable suhordinators. Then for the 
iterated composition defined by Di'^\t) = DloDlo- ■ -oD^^t), the pmf pk{t) = P{N{Di'^\t)) = 
k), satisfies the following equation 

-^Pk{t) = \[Pk{t) - Pk-i{t)]. (3.4) 

Proof. We prove this by induction. The case n = 1 is obtained after putting 5 = 1/ \/2 and 
7 = in (12. 5p . Let k{x,t) denote the density of Dl{t), which can also be obtained by putting 
6 = l/\/2 and 7 = in (II. 6p . For the case n = 2, 



00 poo 



Pk(t) = / / pk{xi)k{xi,X2)k{x2,t)dxidx2. (3.5) 
Jo Jo 



Integrating by parts and then using (13. ip and (13. 2p . we get 

^Pfc(t) = y J Pk{xi)k{xi,X2)-^k{x2,t)dxidx2 



00 poo 



Pkixi)^^k{xi,X2)k{x2, t)dxidx2 



Jo 



00 poo 



JO 



d 

Pk{xi)- — k{xi,X2)k{x2, t)dxidx2 
0x1 



POO POO ^ 

= - / — Pkixi)k{xi,X2)k{x2,t)dXidX2 

Jo Jo '^•^1 

= \[pk(t) -pk-l{t)]. 

The general case also follows in similar way. □ 
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Remark 3.1. (i) The composition of two stable suhordinators is a stable suhordinator. Let 
Di,D2 be two independent stable subordinators with Laplace exponents CiS^\ i = 1,2. Then 

and hence Di{D2{t)) is a stable subordinator of index /3i/32- This implies that the composition 
Ei{E2(t)) is the inverse stable subordinator of index /3i/32- Hence, by induction, we get 

El'\El'\El'\...El''m) = E^'''it). 
(ii) By Remark 2.1, we see that Di"'\t) is a stable subordinator of index 1/2" and 

Note also also that N(D^"'\t)) is also a Levy Process, as N(t) is Poisson process and 
D^'^\t) is a stable subordinator (by Theorem 30.1 in Sato (1999)) and it has Laplace transform 

^(g-s^(DWW)) ^ g-t(log(A(s)))V2" ^ ^^p^_t{X{e'' - 1))V2"), s > 

where fi{s) = E(e~'^^^^'') = exp(A(e~* — 1)). By Theorem 30.1 in Sato (1999), the Fourier 
exponent of N{D^''\t)) can be written as 

E(e^^^(^'"'(*») = exp(-t(A(e^^ - 1))'/'"), z G M. 



Remark 3.2. Let < /3 = ^ < 1 for k,m relatively prime integers, and let D{t) be a stable 
subordinator < 
a solution of 



subordinator of index /3 with E(e *-^(*)) = e In this case, the density f{x,t) of D{t) is 



am pjk 

— fix, t) = {-or— fix, t), x,t>0, (3.7) 
see Lemma 3.1 in DeBlassie (2004)- Hence, the density k*{x,t) of Di^\t) satisfies 

^k*{x,t) = —k*{x,t), x,t>0, (3.8) 

Alternative proof of Proposition [3TT] Using (13. 8 p one can obtain an equation for pk{t) = 
F{N{Di"\t)) = k) as follows: 

Pk{t)= / pk{x)——k*{x,t)dx 



d 







Pk{x)-—k*{x, t)dx 
^ ^ ' ^ (3.9) 

Pk{x)k*{x,t)\^ - I k*{x,t)—pk{x)dx 







= A[pfc(i) -Pk~i{t)], 
which follows by integration by parts and using (13.11) and (13. 2p . 
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Remark 3.3. Consider the standard Cauchy process C{t) with density function 

q{x,t)= / xeR. (3.10) 

tt[x'' + t^j 

The Cauchy process is a symmetric 13 -stable process with index /3 = 1. Its Fourier transform 
(FT) 

r/izs implies 

—q{u,t) = \u\^q{u,t) = -{iufq{u,t). 

Invert the FT to get 

sznce zs i/ie FT ofd^q/dx^. The pmf qk{t) = P{N{\C{t)\) = k), k >0, t > 0, satisfies 

the equation 

(f 2d 

-j-^qkit) = - V)VW - —-^Pk{x)\^=o 

irt ax ]^2) 

= -A^(l - vYqkit) + ^A(l - v)Pk{x)U=o, 

HI 

where \/fk{t) = fk-i{t), the backward shift operator. 

Let f be in the domain of the generator, — A(l — y); of the Poisson process, so that 

-X{l-S7)f{k) = -X{f{k)-f{k-l)). 
Then u{t, k) = Ek{f{N{\C{t)\))) solves 

(f 2d 

— = -A2(1 - \7)\{k,t) - ;^^«(^'2;)U=o 

= -X'il - S7)Mk,t) + -Ml - v)/(^), (2-^^^ 
u{0,k) = f{k). 

A general result that includes Equation f l3.13p as a special case was first proved in Nane l21^ 

We next look at the subordination to inverse stable subordinators. The inverse stable subor- 
dinator E(t) ( of D(t) of index < /3 < 1) has density 

m{x, t) = ^P{E{t) < x) = - P{D{x) < t)) 

= -— r fiu, l)du = (t//3)/(tx-^/^ i)x-'-'/^ 



using the scaling property of the density f{x,t) = {t/(5)f{tx -"^Z^, l)x ^ (see Bertoin 
(1996)). 

It is clear from fl3.14p that the density m{x,t) is in C°°((0, oo) x (0,oo)). Some additional 
properties of m{x,t) are given in the following Lemma whose proof follows from (13.11) . (13. 2p . 
(I3.14p . and by taking Laplace transforms. 

Lemma 3.1 (Lemma 2.1, Hahn et al (2010)). The density m{x,t) of E{t) satisfies 

(a) limt_j._|_o m(x, t) = 6q{x) in the sense of the topology of the space of tempered distribution 

©'(M); 

(b) \im.^^+Qm{x,t) = ff^, t > 0; 

(c) lim^^oo "^(a;, t) = 0, t> 0; 

(d) e~^^m{x, t)dt = s^~^e~^'^'^ . 

Lemma 3.2 (Lemma 2.2, Hahn et al (2010); Theorem 4.1, Meerschaert and SchefHer (2008)). 
The density m{x, t) is the fundamental solution of 

d^mix,t)_ dmix,t) t-^ ^^^^^ 



dtf^ dx r(l - 13) 

in the sense of tempered distributions. 

Next, we consider the iterated composition of Hit) for the case 5 = 1/V2 and 7 = 0. Define 
H'^it) = Hlo ■ ■ ■ oH*(t), where Hl(t), - ■ ■ , H*{t) are n independent copies of H(t). Let qk(t) = 
P{N{H^{t)) = k). By Equation (13. 6p . H''^{t) is the inverse of Di^\t) with Laplace exponent 

Remark 3.4. Let E(t) be the inverse of a stable subordinator of index 1/m then the density 
m{x,t) of E{t) satisfies 

dm{x,t) _ , d"'m{x,t) ^ ^ ^ q. 



dt dx"" 

\x=0 



' - r(^+^/- , t > 0, = 0, 1, 2, ■ ■ ■ , (m - 2 



g^k i^-u -p. (fc+i). 

^^^"^^ (3.16) 



Qm ljj^(^jr^ t) 



0, t > 0; 



lim — ^^^ = 0, t>0,fc = 0, 1,2,--- ,(m-l). 



dx^ 

The PDF (I3.16P was obtained by Keyantuo and Lizama [14^ . 
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Using the last remark we can easily show 



Theorem 3.1. Let E(t) be inverse stable subordinator of index 0</3 = l/m<l, for 
m = 2,3,4, ■■■ . Then the pmf qk{t) = F{N{E{t)) = k) solves the following DDE 

m— 1 



d 

di 



q,{t) = (-A)-(l - vrqkit) + J2 - VYPkiO) 



f~(m-j)/m 



r(i- 



(3.17) 



Proof of Theorem \3.1[ This follows by using integration by parts and fl3.16p above 



d 
It 



d 

Pk{x)—m{x,t)dx 
(-l)™y pk{x)-^m{x,t)dx 

Qj — 1 ffn—j 



™ ai-l 



^-{l+m~j)/m _ 



+ 

x=0 -^0 



gn 



dx' 



-pk{x)m{x, t)dx 



x=o r(i 

^-(l+m-j)/m 

+ 



(1+m-j) ' 



+ 



dx'' 



-pk{x)m{x, t)dx 



m— 1 



,=or(i-M!p) 

^— (1+m— j)/m 







x=0 



r(i 



-A)™(1 - yY'pk{x)m{x,t)dx 



f-im-j)/m 
m 



x=or(i- 



+ (-Ar(i-vrg.w 



(3.18) 



Note that the terms ^-,-i Pfc(3^)U=o might not be zero, in general. For example, PkiO) = 1 for 
k = 0, PkiO) = for > 1. Also, -§^Pkix)\x=o = for A; > 2, and it is equal to A for A; = 1 and 
is equal to —A for A; = 0. □ 

Let f be in the domain of the generator, — A(l — v); of the Poisson process so that 

-A(l-v)/(^) = -A(/(A;)-/(A;-l)). 

Then u{t, k) = Efc(/(A^(E(t)))) solves 



d 



m—l 



^^{t, k) = (-A)-(l - x/ru{t, ^) + E ( (-^)' (1 - V)V(A;) 

j=l \ 

«(0,fc) = f{k),k> 1 
This was first proved by Baeumer et al. (2009) |^ 



f-(.m~j)/ni 
m 



r(i- 



(3.19) 
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Corollary 3.1. The pmf quit) = F{N{H^{t)) = k) solves the following DDE 

-q,{t) = A^"(l - vr<l,{t) + i (-A)^(l - vyp.(O) ] (2^_,) , (3.20) 

Remcirk 3.5. Let E{t) he inverse stable suhordinator of index < ^ < 1. Then the density 
q^[t) = ¥[N{E{t)) = k) solves (see Meerschaert et al. (2011)). 



dt^ 



Qkit) = -A(l - v)(lk{t)- 



(3.21) 



In particular, when 5 = l/\/2 and 7 = 0, the density ?^(t) = P(iV(i7^(t)) = A;) solves 

^1/2" 



-Qm = -K^-V)qm, ri>l. 



(3.22) 



Consider qkit) = P(^N[H"-[E{t))) = k). Since H"{E{t)) is inverse stable suhordinator of 
index /3/2", qk{t) solves 



^0/2" 
dt^' 



■qk{t) = -A(l - V)qk{t), < ^ < 1, n > 1. 



(3.23) 



Arguments similar to the ones above lead to the following result. 
Proposition 3.2. The pmf qk{t) = F{N{H''{E{t)) = k) solves 



dP_ 
dtf^ 



2"-l 



q,{t) = X'"{l-Vrqk{t)+J2 

i=i 



z=0 



■^-/3(2"-,-)/2"^((2"- j)/2^ 

. r(i - (^) 



(3.24) 



where U{-^) = ¥.[D{1)^P] /or < 7 < 1. 
Proof. We can write 

poo poo 

qkit) = ¥{NiH:iEit))) = A;) = / ¥{NiH:ix)) = k)mix, t)dx = / qlix)mix, t)dx. 

Jo Jo 

(3.25) 

Using the fact that the density m(x, t) is the fundamental solution of 

d^mix,t) dmix,t) t~^ 



dt^ 



dx 



r(i-/3) 



Soix), 



(3.26) 
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in the sense of tempered distributions, we get 

d^m{x, t) 



dm{x,t) t~ 
dx 

^ dm{x, t) 
dx 



r(i-/3) 



6o{x) 



dx 



dx - qk{0) 



r(l-/3) 



^dqlix) 
dx 



m(x, t)dx 



+ / ■m{x,t) 
'o 



2"-l 



^-(2"-j)/2" 



(3.27) 



2=0 

oo 



rfi- 



dx 



2"-l 



9^' 



+ E ( ^^^(-) 



=0/ ■>'o 

2"-l 



/3) 

^--(2"-i)/2" 

r(i-^; 



2=0 



^" / m[x,t)qk[x)dx 
Jo 

m{x, t)dx 

(E(t))-(2"-i)/2"- 



E 



r(i 



2" 



We can calculate the terms ^[E{t)-^] for < 7 < 1 as follows: First E{t) = {D{l)/t)-l^ by 
Corollary 3.1 Meerschaert and SchefHer (2004), hence 



00 



Hence, we have 



2"-l 



2=0 



^-/3(2"-i)/2"^((2" -j)/2") 



rfi 



(2"-j)- 



. (3.28) 



□ 



4 Inverse of tempered stable processes as time-changes 

Note that tempered stable processes are obtained by exponential tempering in the distribution 
of stable processes, see Rosinski [22] for more details on tempering stable processes. The ad- 
vantage of tempered stable process over stable process is that they are also infinitely divisible. 



13 



and they have moments of all order. Let f{x,t), < /3 < 1 denotes the density of a stable 
process (stable subordinator) with LT 

POO 

/ e-'''f{x,t)dx = e-''\ (4.1) 
A tempered stable subordinator D^{t) has a density 

^(a:,t) = e-'^^+^''V(x,t), /i > 0. (4.2) 

The Levy measure corresponding to a tempered stable process is given by (see e.g. Cont and 
Tankov (2004, p. 115)) 



ce 



-fix 



^D,{x) = c> 0, X > 0, (4.3) 



oo 



which implies / -K£)^(x)dx = oo and hence using Theorem 21.3 of Sato (1999), the sample 
Jo 

paths of D^{t) are strictly increasing since jumping times are dense in (0, oo). Further the LT 

Us, t) = / e-'^fix, t)dx = e"*«^+^)'-'^'). (4.4) 







Proposition 4.1. The density function f^{x, t) of the tempered stable process of D^{t) satisfies 
the following PDF. For /3 = — , m > 2 



m 

m 



.7 = 1 ^ ^ 



d 

,fiyx,t) = —f^{x,t). (4.5) 

Proof. We prove this result by induction. From (14. 4 p we have 

f^{s,t) = e-"^^^+^^'-^'\ (4.6) 



For m = 2, 



and 



p,is, t) = -((. + i^y/' - /i^/^)/,(s, t) (4.7) 



i2 



^Us,t) = {{s + ^^f'-^,'/'ff,{s,t). (4. 
Using (1121), (B and using the fact that /^(0,t) = (see ([SJD and (jOJ), we get 

^ -2/ii/2|') /,(.,t) = ./,(., t) - /,(0,t). 



Inverting the LT to get 

V/'^)/.(^,t) = ^/,(.,t). (4.9) 
14 



g2 d\ ~ d ~ 



Simililarily, for m = 3 

The result also follows similarily for a general m. □ 
Remark 4.1. Let rk{t) = P{N{Df^{t)) = k). Then the pmf rk{t) satisfies 

r^'''^'"^^"'^^^^ = A(i - vK.(t). (4.11) 

j=i V J / 

Using the steps as the proof of Proposition 4.1., we can show that the density function m^(x, t) 
of the hitting time process E^{t) of Dfj_(t) satisfies, for /? = ^ and m > 2, 

'p^(-iy(^y(^-^/-^)-^rn,ix,t) = ^^m,{x,t)+5oit)m,{x,0). (4.12) 
When /i = 0, (KT2\i reduces to 

Qm Q 

(-l)'"^mo(x, t) = g-^rno{x, t) + 6o{t)mo{x, 0) (4.13) 

Let Efj,{t) denote the inverse of the tempered stable subordinator Dn{t). By Theorem 3.1 in 
Meerschaert and Schefiler (2008) we have the density of E^{t) as 



mf,(x,t)= nD^it-y,oo)ff,{y,x)dy. 
Jo 



Since the density fij.{y,x) of D^{x) is infinitely different iable, m^{x,t) is also infinitely differ- 
entiable. 

Then we have the following result. 
Proposition 4.2. The pmf f kit) = P{N{Ef^{t)) = k) satisfies the following PDE: 

m j / s 

i+fc/"^ \ ,,(l-i/m) 



7 = 1 V"' / 7 = 1 k=l 



j=i fc=i ' (4.14) 

X ^-^ -Pk{x)-^. -m^(x,t) -5o(t)rfc(0). 

oxk — 1 c/xj — /c lx=o 
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Proof. We have 



d 

Pkix)—m^{x,t)dx 



m 



dt 



J 







■rn J 



^(_i).(";)^(i-.7H 



r 5 S-'" 



fc=i 



a;=0 



d^ 



+ {-^y I -^Pk{x)m^{x,t)dx -5o{t)rk{0) 



j=i k=i 



Qk- 



i-k 



x=0 



Hence the result. 



□ 
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5 Notation 



variable 


density 


G{t) 


g{xj) 


S = l/\/2,7 = 


k{x, t) 


H{t) 


h{x, t) 


5 = 1/A7 = 


l{x,t) 


/3 - stable D{t) 


fM 


13/2- stabler (t) 


p(x,t) 


1/2 — stable n* 


k*{x,t) 


m 


m{x, t) 


6 = 1/^2,7 = 


l{x, t) 


1/2 — inverse stable n* 


l*{x,t) 


Caucy C{t) 


q{x,t) 


TSD^it) 


Mx,t) 


Hitting time E^{t) 


m^{x,t) 
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